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Abstract. Substances such as the "telephone number compound" Sri4Cu24 04i are intrinsically hole-doped. 
The involved interplay of spin and charge dynamics is a challenge for theory. In this article we propose to 
describe hole-doped Heisenberg spin rings by means of complete numerical diagonalization of a Heisenberg 
Hamiltonian that depends parametrically on hole positions and includes the screened Coulomb interaction 
among the holes. It is demonstrated that key observables like magnetic susceptibility, specific heat, and 
inelastic neutron scattering cross section depend sensitively on the dielectric constant of the screened 
Coulomb potential. 

PACS. 75.10.Pq Spin chain models - 75.40.Mg Numerical simulation studies 



1 Introduction and model 



Substances hosting spin and charge degrees of freedom 
^ exhibit a large variety of phenomena like magnetic and 

■ charge ordering, metallic conductivity and superconduc- 
I tivity PE]. The "telephone number compound", 

. Sri4Cu2404i, contains two magnetic one-dimensional struc- 
' tures, chains and ladders. The stoichiometric formula sug- 
, gests 6 holes per formula unit. We will assume that for the 

■ undoped compound all holes are located in the chain sub- 
' system (i. e. 60 % holes), although this is experimentally 
, under discussion since x-ray absorption (XAS) measure- 
' ments suggest that at room temperature some holes are 
I located in the ladder subsystem [H], whereas it is neces- 

■ sary to assume that all holes are in the chain subsystem 
' in order to explain neutron scattering data HI. 

At low temperatures the ladder subsystem is magnet- 
' ically inactive due to a large spin gap [5|. The remaining 
, dynamics of the hole-doped chain system is still interest- 
ing as well as complicated enough to constitute a challenge 
] for theoretical investigations. Especially the evaluation of 
1 thermodynamic quantities both as function of tempera- 
ture and magnetic field is prohibitively complicated even 
for moderate system sizes. Therefore, mostly approximate 
descriptions in terms of classical spin dynamics ||6i|7| , spin- 
dimer models |51IMTn| or spin- wave analysis ^Jj have been 
applied. Calculations based on the Hubbard model aim at 
ground-state correlations at low hole doping |12j . 

A fundamental question in this context is how the 
charge order in the Cu02 chains of substances such as 
Sri4Cu2404i is established. These chains seem to be a 
sequence of rather perfect antiferromagnetically coupled 
spin dimers separated by holes, see Fig.^ Any proposed 
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theoretical model should also be able to describe excita- 
tions involving hole motion which is crucial because inter- 
esting physical properties of these compounds result from 
a competition of charge mobility and magnetic interac- 
tions |13II14II15II16| . One possible explanation is that the 
formation of dimers is generated by structural modula- 
tions of the material via a strong variation of the on-site 
orbital energies |17II18II19II20| . 

In this article we investigate how a screened electro- 
static hole-hole repulsion along the chain would express 
itself in thermodynamic quantities. It will turn out that a 
rather strong Coulomb repulsion is needed in order to re- 
produce the experimental magnetization. This is in accord 
with e.g. Ref. [2] or with Refs. j22EBl ■ where a dielectric 
constant of 3.3 is found to be realistic. 

In order to be able to evaluate thermodynamic quan- 
tities we propose to describe hole-doped spin rings with a 
Heisenberg Hamiltonian that depends parametrically on 
hole positions. This ansatz is similar to a simple Born- 
Oppenheimer description where the electronic Hamilto- 
nian (here spin Hamiltonian) depends parametrically on 
the positions of the classical nuclei (here hole positions) . 
Each configuration c of holes and spins defines a Hilbert 
space which is orthogonal to all Hilbert spaces arising from 
different configurations. The Hamilton operator H{c) of a 
certain configuration c is of Heisenberg type and depends 
parametrically on the actual configuration c, i. e. 

H = J2 {Hie) + Vic)) (1) 

C 

Hie) = -Y, Juvic) siu) ■ siv) . (2) 

u,v 

Juv{c) are the respective exchange parameters which de- 
pend on the configuration of holes. J < describes anti- 
ferromagnetic coupling, J > ferromagnetic coupling. 
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Fig. 1. L.h.s.: Ground-state hole configuration for 20 sites and 
60% holes. This configuration is also called dimer configuration 
since it consists of weekly interacting antiferromagnetic dimers. 
One dimer is highlighted (d). The hole-spin exchange processes 
1 and 2 lead to energetically low- lying configurations. R.h.s.: 
Exchange parameters used in this article: J = —67 K, Jy = 
7 K, and Jnn = 25 K. 

Figure^shows on the l.h.s. as an example the ground- 
state hole configuration of Sri4Cu2404i pilllll^ which 
is a sequence of spin-hole-spin dimers separated by two 
holes. Energetically excited configurations arise if holes 
are moved to other sites as depicted exemplarily by the 
exchange processes 1 and 2. The r.h.s. of Fig. ^ illustrates 
how the exchange parameters depend on the actual hole 
configuration. In this work three different exchange pa- 
rameters are employed. 

A key ingredient of the proposed model is the inclu- 
sion of the electrostatic interaction between holes which 
is modeled by a screened Coulomb potential 

where tq — 2.75 A is the distance between nearest neigh- 
bor sites on the ring. The dielectric constant is con- 
sidered as the only free parameter in the present inves- 
tigation. Several attempts have been undertaken to esti- 
mate the dielectric constant which yielded values for up 
to 30 P^25ll2Hj . In related projects where the exchange 
interaction of chain systems in cuprates is derived from 
hopping matrix elements between different orbitals using 
a Madelung potential the dielectric constant is found to 
be Er 3.3 EUSl- 



of hole configurations is already too big to be considered 
completely. Therefore, only the ground-state configuration 
and low-lying excitations with their respective degenera- 
cies are taken into account. This is sufficient since hole 
configurations which deviate considerably from the ground 
state configuration possess very high excitation energies. 

Correlated electrons are usually modeled with the Hub- 
bard model J therefore looking at Eqs. and |j2Jl one 
might be tempted to ask: Where is the kinetic energy of 
the holes? The Hamiltonian of the Hubbard model pTifeHl 

transforms at large U and half filling into a Heisenberg 
Hamiltonian (O with Jij = —4tfj/U. Therefore, if working 
close to half filling, it is legitimate to say that the kinetic 
energy is absorbed into the exchange coupling. Neverthe- 
less, the remaining hole motion is treated classically, i. e. 
superpositions of hole states are not taken into account. 
For the actual compound which is far from superconduc- 
tivity this assumption of well-localized holes seems to be 
appropriate. 

What is properly taken into account is the screened 
Coulomb interaction between holes. But, if so: Wouldn't 
it be sufficient to consider nearest-neighbor Coulomb re- 
pulsion only? Although this is a common method wc find 
that a simple nearest-neighbor repulsion results in unphys- 
ical ground states. It is experimentally verified by means 
of low-temperature susceptibility |l(Jj , neutron scattering 
1^1 II as well as thermal expansion measurements 9 , that 
the highly symmetric dimer configuration, compare Fig.^ 
constitutes the ground state of Sri4Cu2404i. Using only 
nearest-neighbor Coulomb repulsion yields an alternating 
sequence of spins and holes with the remaining 10 % holes 
assembling as a big cluster irrespective how strong the 
repulsion is. The reason is that this strange configuration 
has the same number of nearest hole-hole neighbors as the 
dimer configuration. Even the inclusion of a next-nearest 
neighbor Coulomb repulsion does not improve the situa- 
tion, the Coulomb interaction is still proportional to the 
number of sites and may be overcome by the antiferro- 
magnetic binding J. 



2 Discussion of the model 

The aim of the proposed model is to evaluate the complete 
spectrum for reasonably large system sizes and thus to 
be able to investigate thermodynamic quantities both as 
function of temperature and field. The spectrum does not 
only consist of those levels belonging to the ground-state 
hole distribution, compare Fig.^ but also of all levels aris- 
ing from all other hole configurations. For small systems 
all such configurations can be generated and the related 
spin Hamiltonians (0) can be diagonalized completely. For 
8 spins and 12 holes for instance this amounts to 6310 
distinct hole configurations and tiny Hilbert spaces of di- 
mension 256. For 16 spins and 24 holes the total number 



3 Results 

At temperatures and energies below 200 K the behavior of 
the chain subsystem in Sri4Cu2404i is usually discussed in 
terms of weakly interacting dimers sometimes augmented 
by weak interchain interactions, see e. g. |Hllll','^T0'. Such 
a picture, although rather successful, does not allow to 
discuss the influence of mobile holes on thermodynamic 
observables. It is clear that configurations like 1 and 2 in 
Fig. n will contribute to thermal averages, but how? 

Figure [5] shows the low-energy part of the spectrum 
of a chain of A^tot — 40 sites with A^s = 16 spins and 
A^h = 24 holes for three choices of the dielectric constant 
Er. If Ef. = 1 the spectrum up to several hundreds of Kelvin 
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Fig. 2. Low-energy part of the spectrum of a chain of A'tot = 40 
sites with A's = 16 spins and Nh ~ 24 holes for three choices 
of the dielectric constant e^: bars - er — 1, bars together with 
x-symbols - tr = 5, bars and crosses - er — 10. The arrow 
marks the singlet-triplet transition employed in dimer models. 



is solely given by the levels of the ground-state dimer con- 
figuration (bars in Fig.|2). With increasing Cr the Coulomb 
repulsion decreases and so does the excitation energy of 
magnetic levels belonging to hole configurations where one 
or two holes are moved. As an example the levels resulting 
from such configurations are given as x-symbols (e^ = 5) 
and crosses (e^ — 10) in Fig.|21 It is clear that besides the 
singlet-triplet transition (arrow in Fig. [SJ, which is the 
main ingredient of the dimer model, transitions to states 
involving spin-holes exchange processes will contribute to 
thermodynamic observables like the inelastic neutron scat- 
tering cross section. In the following we discuss the influ- 
ence on three basic observables. 
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Fig. 3. Zero- field magnetic susceptibility xoiT): solid curve - 
e,. = 1, dashed curve - tr = 5, dashed-dotted curve - er = 10. 
For a comparison experimental data, taken at B = 1 T, are pro- 
vided by x-symbols 1101 . The corrected data, given by crosses, 
take also into account that due to impurities the number of 
dimers is less than theoretically possible |10) . 



effect is at all observable since the responsible levels are 
at excitation energies well above that temperature range, 
compare the spectrum in Fig. |2] 

It turns out that the high degeneracy of excited hole 
configurations is the reason for the influence of the hole 
dynamics even at low temperatures. Looking again at the 
hole configuration shown in Fig. ^ one notices that the 
spin-hole exchange processes can happen at very differ- 
ent places leading to a large geometric degeneracy. This 
degeneracy can overcompensate a small Boltzmann factor 
and thus expresses itself in a high thermal weight. 

Together with the theoretical results Fig. 01 shows the 
experimentally determined magnetization which was mea- 
sured along the c-axis of the material at a magnetic field 
of B = 1 T jJU). The x-symbols depict the uncorrected 
values whereas the crosses represent the corrected values. 
The correction includes a subtraction of impurities (free 
spins s — 1/2) as well as a rescaling because the number 
of dimers on the chain is less than theoretically possible. 
The almost perfect coincidence with the theoretical re- 
sult for er — I suggests that the hole-hole repulsion is 
rather strong. The uncertainties in the measurement and 
the correction procedure leave some freedom for the ac- 
tual value, but it is clear that should not be bigger than 
about three. This implies that energy levels which result 
from other than the dimer configuration are well above 
the triplet excitation, compare Fig. [5] Although this might 
seem to be unrealistic one has to keep in mind that any 
theoretical model must explain why the experimental sus- 
ceptibility practically coincides with that of free dimers. 
This is only possible if other excitations are well separated 
from the triplet excitation. 




Fig. 4. Specific heat at _B = 0: er = 1 - solid curve, er — 2 - 
dotted curve, £r = 3 - thin curve, = 5 - dashed curve, and 
Er = 10 - dashed-dotted curve. 



Figure 131 presents the results for the magnetic suscep- 
tibility x{T,B = 0) = Xo(r) at vanishing magnetic field 
B — 0. The solid curve shows the theoretical susceptibility 
for Er = I, the dashed curve for = 5, and the dashed- 
dotted curve for e^- = 10. One realizes that with increasing 
Cr, i. e. with stronger screening of the Coulomb interaction, 
the susceptibility increases at intermediate temperatures 
and that the maximum shifts to lower temperatures. Al- 
though being rather moderate it is astonishing that the 



But even if higher-lying energy levels are well above 
the triplet excitation, due to their vast degeneracy they 
can substantially contribute to the specific heat, which is 
shown in Fig. 0] In order to demonstrate how the ther- 
mal weight of the excited hole configurations grows sev- 
eral cases are shown. The solid curve, which is the lowest 
among all curves, again depicts the result for er — I. With 
increasing dielectric constants thermal weight is shifted 
from higher temperatures to lower ones. This can very 
clearly be seen for the shown sequence: er ~ 2 (dotted 
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curve), er — 3 (thin curve), = 5 (dashed curve), and 
€r = 10 (dashed-dotted curve). This result shows that for 
dielectric constants of the order of ~ 5 about one third 
of the specific heat at its maximum is due to states in- 
volving hole motion. 

Inelastic neutron scattering is a valuable tool to mea- 
sure the magnetic excitation spectrum of substances like 
Sri4Cu2404i, see e. g. pil3(J) . The transition from the sin- 
glet ground state to the first excited triplet state, arrow 
in Fig. 121 has been measured with high accuracy and used 
to determine the exchange constants, especially J (Fig.^ 
r.h.s.). For recent results have a look at . 




E (K) 

Fig. 5. Rough sketch of the lowest transitions observable with 
inelastic neutron scattering: er = 1 - solid curve, er = 5 - 
dashed curve, and = 10 - dashed-dotted curve. The arrow 
marks the singlet-triplet transition at about 135 K. 

In addition to this fundamental transition of the dimer 
configuration, transitions to states with spin-hole exchange 
should be detectable, too. FigureElshows as a rough sketch 
where such transitions could be expected. The singlet- 
triplet transition at about 135 K - arrow in Fig. |5l - is 
clearly seen for each dielectric constant. But with increas- 
ing er transitions to states with one or two holes moved 
become accessible. The dashed curve shows schematically 
the lowest transitions for = 5, the dashed-dotted curve 
the lowest transitions if = 10. Although the singlet- 
triplet transition might have the largest matrix element, 
again the huge geometric degeneracy could help to make 
the other transitions visible. A possible drawback is, nev- 
ertheless, due to the fact that states involving spin-hole ex- 
change break the translational symmetry, compare Fig. ^ 
The momentum dependency, therefore, might be fuzzy. 

4 Summary and Outlook 

For the compound under investigation the model is ca- 
pable to address much more questions. Since it is not 
clear whether the hole content of the chain subsystem in 
Sri4Cu2404i is really 60 %, one can study the influence 
of a reduced number of holes on magnetic observables. 

Figure shows as an example how a reduction by just 
one hole will influence the magnetic susceptibility. The 
perfect symmetry of successive spin- hole-spin dimers sepa- 
rated by two holes is destroyed and a low- lying triplet com- 



20 




Er=6, :=(-67, 7, 


25) K, B=0 


15 


1 






10 








0,5 




/ N," 1 6 , 


Nh = Z4 






/ N , = 1 6 ,, 


N,= 23 


00 













50 100 


150 200 






T (K) 





Fig. 6. Zero-field magnetic susceptibility Xo{T) for e,, = 5: 
solid curve - A's = 16 and A'h = 24, dashed curve - A'^g = 16 
and iVh = 23. 

petes with the former singlet ground state which leads to a 
low-temperature divergence of the susceptibility. It might 
very well be that a part of the experimentally observed 
low-temperature divergence of the susceptibility |l()j is due 
to the reduced hole doping of the chain. 

The absence of perfect symmetry also leads to an in- 
creased mobility of the holes. At the imperfections neigh- 
boring holes can be moved without altering the Coulomb 
energy much. Therefore, the excitation energy for config- 
urations where a hole is moved at imperfections will be 
rather low. 

Summarizing, the main advantage of the proposed ef- 
fective spin-hole Hamiltonian is that it allows to evaluate 
thermodynamic observables both as function of tempera- 
ture and magnetic field for reasonably large systems. Us- 
ing this model it could be shown how a hole-hole Coulomb 
repulsion along the chain would express itself in thermo- 
dynamic observables. The comparison with experimental 
magnetization data suggests that the screening is weak. 
The actual choice of the exchange parameters, compare 
Fig. m does not influence the general conclusions about 
the effects of the Coulomb repulsion between holes. It is 
not yet clear whether a weekly screened Coulomb repul- 
sion is realistic PI] or whether a combination of hole-hole 
Coulomb repulsion and modulation of on-site energies |17l 
J,8, 19 ,20^ has to be used. An experimental determination 
of other excited states than the triplet excitation of the 
dimers would be very helpful in this respect. 
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